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Many independent high resolution simulations have indicated that the standard collisionless cold 
dark matter model does not reproduce the structure of observed present day galaxies well. Several 
possible solutions in the form of modifications to the physics of the dark matter particles have 
been proposed. One of the most promising is warm dark matter (WDM), particles with significant 
thermal motion in the early universe. It is usually assumed that such particles are relativistically 
decoupled particles with a mass of approximately 1 keV. However, here we have investigated the 
possibility that much more massive particles with highly non-thermal spectra could make up warm 
dark matter. Several possible production mechanisms are reviewed and the only one found to be 
viable is that the WDM is produced by the non-relativistic decay of some massive species in the 
early universe. Such very massive warm dark matter could possibly be detected in direct detection 
experiments, as opposed to standard thermal warm dark matter. 

PACS numbers: 95.35. +d, 98.65. Dx, 14.80.-j 
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I. INTRODUCTION 

The standard Big Bang model has been very successful 
in explaining the large scale structure in the universe. An 
essential feature in the model is dark matter. By far the 
most successful candidate for dark matter is cold dark 
matter, collisionless particles which are so massive that 
they are very non-relativistic during the entire structure 
formation history . 

However, in the past few years high resolution N-body 
simulations have shown that standard collisionless CDM 
apparently fails to reproduce observations on galactic or 
smaller scales. Observations indicate that the halo of a 
galaxy like the milky way has almost an order of mag- 
nitude fewer small satellite halos than is found in CDM 
simulations 

Also, the central cores of dark matter halos have sin- 
gular density profiles, approaching almost p oc r -1 ' 5 at 
small r [B-Bfl. This is in contrast to observations of low 
surface brightness galaxies, where the dark matter halos 
apparently have almost constant density. 

This seems like a very serious problem for the CDM 
structure formation picture. However, at present it can- 
not be completely excluded that the discrepancy is due 
to the low quality of either simulations or observations. 
For instance very high resolution simulations which in- 
clude baryons have yet to be carried out , and rotation 
curve measurements of LSB galaxies may have greater 
uncertainties than previously estimated 

However, the amount of observational data in conflict 
with the CDM model is quite large and it seems a very 
real possibility that we may have to abandon the sim- 
plest CDM collisionless models. The simplest possibility 
(other possibilities, such as self-interacting CDM have 
been proposed [Q) is that power is suppressed on small 
scales, either because of the initial power spectrum from 
inflation 111311 or because of thermal motion of the dark 



matter particles. The latter possibility corresponds to 
warm |[4|,[l5) or hot dark matter and has received a lot of 
attention recently. WDM seems able to explain many of 
the problems which CDM suffers from, notably the sub- 
structure and the angular momentum problems jl6] E9[. 
However, it remains to be seen whether WDM can also 
prevent singular cores from forming fll8| , p0| . 

Because WDM seems a quite promising candidate for 
dark matter it is definitely worthwhile to look closer at 
the possible particle physics mechanisms for producing 
it. The simplest possibility is that the warm dark mat- 
ter is a relativistically decoupled species with a thermal 
distribution function. This possibility has been reviewed 
a number of times and we shall not discuss it further 
@,|5||l). Rather we will investigate the intriguing pos- 
sibility that warm dark matter could have been produced 
non-thermally. It is possible that dark matter particles 
with masses much higher than thermal warm dark matter 
could have been produced with enough thermal energy to 
make up warm dark matter. We shall discuss three dis- 
tinct possibilities for this, namely production by decay 
of a heavy species [fl8|| , production by annihilation of a 
heavy species, and finally self heating by number chang- 
ing self interactions of the warm dark matter particles 
p2|- p4] , |l9|] . We then go on to discuss any possible sig- 
natures that can distinguish massive warm dark matter 
from ordinary warm dark matter. 

Note that there is another possibility for non-thermal 
warm dark matter, namely that the WDM particles are 
very light, but have been produced with very low average 
momentum. This could happen if the WDM particles are 
sterile neutrinos that have been produced via oscillations 
with an active species pq| or byrelativistic decays of 



massive fermions into bosons 26 2 
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II. PRODUCTION MECHANISMS FOR 
MASSIVE WARM DARK MATTER 

A non-thermal warm dark matter candidate could have 
been produced in several different ways via interactions 
with other particles in the plasma. We here review three 
distinct possibilities: self-interactions, decays and pair- 
annihilations. In the following, cj> is used to denote either 
the WDM field or the WDM particles, irrespective of 
whether they are bosons or fermions. 



A. Production by self interactions 



Carlson et al. |22 23 suggested an interesting way to 
heat cold dark matter relative to standard CDM, namely 
by number changing self-interactions. This was pro- 
posed as a means of reconciling the Sl m = 1 CDM model 
with large scale structure observations. However, subse- 
quently it was shown by de Laix, Scherrer and Schaefer 
[pifl that this model produced a very poor fit to observa- 
tions. 

At present the consensus is that the correct model for 
large scale structure is a flat ACDM model |jl|, and so 
there is no need to modify the physics of CDM to change 
the formation of large scale structure. Self-interacting 
dark matter could instead be invoked as a possible means 
of explaining the observed small scale structure. CDM 
with a large self-scattering cross section has been pro- 
posed Ip^l , and subsequently warm dark matter with 
self-scatterings |2q| . However, if the self-interactions are 
sufficiently strong there are very likely number chang- 
ing reactions as well. This model is almost equivalent 
to that originally proposed by Carlson et al. ^2j,^3|. If 
dark matter has such strong self-interactions that num- 
ber changing reactions occur in equilibrium, the effective 
temperature of the species will drop much slower than for 
other species. In case of full thermal equilibrium main- 
tained by self-interactions, the distribution function is 
(assuming Boltzmann statistics) / = exp(— E/T). Also, 
in thermal equilibrium the total entropy within a comov- 
ing volume is conserved. In this case, the effective tem- 
perature of the distribution is 



- = (31og(a) 
m 



K)~ 



(1) 



for T <C m. a is the scale factor and K is a constant. 
This shows that the effective temperature of the species 
only drops logarithmically. Note that this behaviour 
is in stark contrast with that of the other species for 
which T oc a -1 . Thus, the particle distribution is heated 
to higher temperature than the surrounding medium by 
transforming particles into kinetic energy. This leads to a 
much larger free streaming length for such particles than 
for thermally decoupled particles with the same mass. 
The observational bounds on the self-scattering cross sec- 
tion can be translated into a rough bound on the num- 



ber changing reactions, depending on the specific type of 
WDM particle. 

Pseudoscalars — A likely candidate in this category is 
a particle like the majoron (although there should be no 
coupling to standard model fields) with a simple atfi self- 
interaction term. 4^2 number changing interactions 
take place via a (f> 6 term in the effective Lagrangian (^|| . 
Based on purely dimensional arguments the rate for this 
reaction should be 1241 



n a /m 



(2) 



where n denotes the number density. In all reasonable 
scenarios we should expect that the dimensionless cou- 
pling constant, a < 1. However, we can also put another 
constraint on a from simple elastic scattering. The rate 
for elastic scattering of these particles is 



2 1^,2 

nva /m , 



(3) 



where v is their relative velocity. It has recently been 
shown in numerical simulations that dark matter with a 
large elastic cross section is not a good candidate because 
galaxies have very singular halos and are too spherical 
fll2| ] . A safe upper bound on the 2-body elastic scattering 
cross section is [n2| 



a < 10 22 mGevcm 2 , 



yielding an upper limit on a 
a 2 < 2.6 x 10" 



22 „ 



(4) 



(5) 



The quantity of interest here is the ratio of the num- 
ber changing interaction rate and the Hubble parameter, 
T/H. For pseudo-scalars this can be written as 

I < 1.1 x ltf\n^fa*m:?Tl Y g- l l\ (6) 

using the bound on a. f2^ is the present-day density 
contribution of the WDM and g* is the number of rela- 
tivists degrees of freedom present. Freeze-out occurs at 
T/H ~ 1, corresponding to the bound 



5/7 < 



"eV 



0.0017(^ 2 ) 3 / 7 T frcczc 



cV- 



(7) 



As soon as the number changing reactions freeze out the 
particles become non-relativistic. Therefore, if this sce- 
nario is to work, we must demand that Tf roozc < 100 — 200 
eV, so that the WDM has the correct amount of thermal 
motion. For pseudo-scalars this translates into 



m cV 



< 0.22(^ 2 ) 3 / 5 « T frcGZCi 



eV- 



(8) 



The conclusion is that massive WDM cannot be produced 
in this fashion unless special circumstances prevail. One 
way out is to tinker with the 4^2 reaction rate, increas- 
ing it by a large factor over the naive estimate. However, 
it seems safe to assume that this type of heating is an 
unlikely production mechanism for WDM. 
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Scalars — If the particles are scalars there could be a 
«-» 2 number changing term of the form T ~ n a 3 /m 5 
||. In this case Eq. (||) changes into 



m cV < 4.1(ft /i 2 ) 2 / 5 



(9) 



Again, the conclusion is that self-heating cannot be the 
production mechanism. 

Fermions — For fermions the situation is quite differ- 
ent. The leading number changing reaction is 2/ + 2/ — » 
//. The fcrmion-fermion interaction necessitates some 
new force carrier, 77. If it is a boson with mass smaller 
than the fermion then the reaction // — > 2r] will dom- 
inate completely and the fermion will quickly disappear 
by annihilation. Thus, the force carrier would have to 
be more massive than /. But in that case the reac- 
tion 2/4-2/ — > // will be suppressed by high orders 
of T/m n <C 1. Thus we can expect the number chang- 
ing rate for fermions to be much smaller than for pseu- 
doscalars and far too small to be of any practical interest. 

Thus, no matter what nature the dark matter particles 
have, it is very unlikely that self-heating can increase the 
thermal energy sufficiently that WDM results. 



B. Production by decay of a massive species 

Another, and perhaps more obvious way is to pro- 
duce the WDM particles,<£, via decays of strongly non- 
relativistic particles, X [Q. We can write the present 
day thermal velocity of the warm dark matter as 



m<p to^ \g*,d 



1/3 



(10) 



AN„ = 0.1850*^ 



0.4km s 



(13) 



For = 0.3 and h = 0.65 this yields AJV„ ~ 0.023, too 
small a perturbation to be detected. Even with high pre- 
cision measurements of the primordial abundances this 
value seems out of reach If the decay is after BBN, 
the contribution of X during BBN is entirely negligible. 

Next, one can ask whether this is detectable in the 
CMBR anisotropy. At recombination the ^-particles are 
already strongly non-relativistic. Their contribution to 
the radiation energy density can be estimated as 



P<f>,R 



Iv^h 2 ( 



(14) 



Again, we can parameterize this, in units of neutrino 
species at recombination, to be 



AJV„ 



1.7 x 10- 



0.4km s 



T=T d , 



At recombination this gives AN V ~ 2 x 10~ 5 , which is 
too small to be detected, even with the upcoming high- 
precision experiments MAP and PLANCK 30 1. 

Note that thermally decoupled WDM gives the same 
contribution to the cosmic radiation density as does the 
decay-produced WDM. Thus, this effect would in no case 
allow us to distinguish the two different types of dark 
matter (unless the decay is after BBN, in which case the 
decay-produced WDM would not contribute to the radi- 
ation during BBN). 



C. Production by annihilation of a massive species 



We can use this relation to find how much energy density 
was in the X-particles when they decayed 



PX,d = Pc^cf 

m 4> \9*.oJ \-L~f,o 



(11) 



where 0,$ is the present day contribution to from the 
warm dark matter. From these relations we can calculate 
the ratio px / Pr at decay 



— = 0.025^/i 2 
PR 



0.4km s 



\ 1/3 

g*A (12) 

9*,oJ 



From this we can conclude that the A-particles never 
dominated the energy density of the universe prior to 
decay. However, the radiation produced at decay adds a 
significant contribution to pa which could be detectable. 

If the decay took place prior to BBN (z ~ 10 12 ) the 
<p particles were highly relativistic during BBN and their 
contribution to the radiation energy density is (in units 
of the energy density of a standard massless neutrino 
species) 



The final production mechanism we shall discuss is via 
the pair-annihilation of some massive fermion species. 
Here, the ^-particles are produced via XX — * 20. <f> 
can be either a boson or a fermion. The prime exam- 
ple would be a heavy neutrino species annihilating into 
massless fermions. Annihilations will in general not pro- 
duce 4> with much greater than thermal energies. The 
reason is that T/H is always an increasing function of T 
for non-relativistic particles, so that very few high energy 
annihilation products are produced (unless one invokes a 
specific physical mechanism to prevent annihilations at 
high temperature |pl[ ). 

In any case, if warm dark matter is produced by the an- 
nihilation of a massive species, the annihilation rate must 
be high enough that practically all heavy particles disap- 
pear before the annihilation freezes out. Otherwise the 
remaining heavy particles will quickly begin to dominate 
the energy density, resulting in a standard CDM scenario. 
But in the limit of strong interactions, the specifics of the 
annihilation reaction does not matter, the decay proceeds 
in equilibrium. The final energy density of the annihila- 
tion products can then be found by solving the simple 
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differential equations (assuming 2 <-> 2 annihilation and 
Boltzmann statistics of the involved species) 



dnxoT dnx 
~dT~ ~ ~dT 4 

dp-TOT _ dpx + 
dt 



dn^ 
~~dt 



-3H TlTOT 



(16) 



-3ff(pTOT+P TO T). (17) 



fx 



ex.p(-{^ml t<p +p 2 -fi)/T), 



(18) 



where and T denote the common pseudo-chemical po- 
tential and temperature of the two species. 

Fig. 1 shows the outcome of the annihilation process. 
The average energy of the light particles is slightly higher 
than the expected 3T, namely {E$) = 3.78T. However, 
this clearly shows that any annihilation-product cannot 
act as strongly non-thermal warm dark matter. 

After reviewing these three possible production sce- 
narios for non-thermal WDM we conclude that the only 
plausible possibility is that the WDM particles have been 
produced by the decay of some massive species. 
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FIG. 1. The evolution of the energy density of parent 
(dashed line) and daughter (full line) for annihilation in equi- 
librium, parametrized in units of the initial high temperature 
density for both species (which is the same since they are in 
thermal equilibrium). Note that the density is multiplied with 
a 4 so that it is constant for a decoupled relativistic species. 



III. COUPLING TO STANDARD MODEL FIELDS 

Warm dark matter could in principle be coupled to 
standard model fields. Thermal dark matter is much 
lighter than the Z-mass and therefore any coupling to the 
standard model fields would likely show up as a branch- 
ing to 4> in the Z-decay. However, thermal WDM cannot 



be detected in standard dark matter search experiments 
that rely on nuclear recoil effects because the mass is 
much below detection threshold That would not 

need be the case for massive warm dark matter candi- 
dates. In principle the WDM particles could scatter on 
nuclei and be observed in direct detection experiments. 
However, the WDM should never come into equilibrium 
with the standard model particles in the early universe. 
Otherwise a purely thermal distribution results, and the 
end product is cold dark matter. This puts a strict up- 
per limit on the coupling between <p and standard model 
fields. 

As an example we calculate the energy equilibration 
rate for a Dirac neutrino-like particle, </>, in the early 
universe. The squared matrix element for (f> scattering 
on a massless standard model fcrmion is roughly 



(19) 



for energy transfers below mz- The dimensionless pa- 
rameter i] denotes the effective coupling strength of <fi 
to Z. From the above matrix element we estimate the 
scattering rate for a </>-particle to be 



r. = n 



r)G%(Et)(E f )Tj. 



(20) 



The average energy of <fi is much bigger than that of /. 
Thus, on average the momentum transfer in each scatter- 
ing is equivalent to (Ef). The energy equilibration rate 
can then be estimated as 



r E ~r)G%(E f ) 2 Tf~r)G 2 F T. 



L f- 



(21) 



If the WDM particles are to stay out of thermal equilib- 
rium then 



F 
H 



< 1 



(22) 



must be fulfilled at all times after <f> is produced. This 
above equation can be transformed into a bound on r\ 



rj < 2.8g*' Tj 



l/2 T -3 



McV 



(23) 



Thus, the bound strengthens with increasing tempera- 
ture and is strongest at the decay temperature of the 
parent. 

A direct detection experiment would typically use 
WIMP-nucleon elastic scattering recoil. The elastic scat- 
tering cross section is given approximately by p2| 



cr ~ 1]G 2 F fl 2 , 



(24) 



where /x = (misrm^)/(mN + m^,) is the reduced mass of 
the system. 

The present generation of experiments have reached a 
limit of cr ~ I pb for masses of 4> in the GeV range 33 1 . In 
that case p/rriN ~ I and the detection limit corresponds 
roughly to 



T) ~ 20. 



(25) 
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Using Eq. ( |23| ) we find that 

Td.MeV ^5 0.5 



(26) 



if the (/(-particles were to have been detected in the 
present experiments. The next generation of dark mat- 
ter experiments will probably achieve about a factor 10 3 
better sensitivity, yielding T^McV ^5 5. In any case, if 
this type of warm dark matter is to be detected in direct 
detection experiments it must have been produced at the 
epoch of BBN or later, otherwise the upper bound on the 
coupling to standard model fields becomes so tight that 
no direct detection experiment will be able to see it. 



IV. DISCUSSION 

We have discussed various possible means of produc- 
ing warm dark matter, with specific focus on the pos- 
sibility of very heavy warm dark matter. Dark matter 
self-heating was found to be excluded by present obser- 
vational bounds. Also, production by pair-annihilation of 
a massive species cannot produce strongly non-thermal 
dark matter. The only viable production mechanism was 
found to be production by the non-relativistic decay of a 
massive relic. 

From a structure formation point of view there is very 
little difference between massive warm dark matter and 
standard thermally decoupled warm dark matter for the 
same average free streaming length. Interestingly , as op- 
posed to standard warm dark matter it could be possible 
to detect massive warm dark matter in a direct detec- 
tion experiment, because the particle mass can easily be 
in the GeV range or above. However, it was also found 
that such WDM would have to have been produced at 
relatively low temperatures, otherwise the WDM distri- 
bution would have been brought into thermal equilibrium 
by interactions with the standard model fields. 
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